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Abstract. Let A be a self-adjoint operator on a separable Hilbert space S). As- 
sume that the spectrum of A consists of two disjoint components ctq and C7i such 
that the set <Jo lies in a finite gap of the set CTi. Let V be a bounded self-adjoint 
operator on S) off-diagonal with respect to the partition spec(A) = cTq U (Ti. It is 
known that if < \/2d, where d = dist(cTo! C^l )> then the perturbation V does 
not close the gaps between Oq and tTi and the spectrum of the perturbed opera- 
tor L = A + V consists of two isolated components coq and ©i originating from 
Go and CTi , respectively. Furthermore, it is known that if V satisfies the stronger 
bound ||V|| < d then for the difference of the spectral projections Eyi((To) and 
Ei{coq) of a and L associated with the spectral sets cTq and coq, respectively, the 
following sharp norm estimate holds: 

/ \\V\ 

l|EA(cfo)-EL(cuo)ll <sm( arctan — 

In the present work we prove that this estimate remains valid and sharp also for 
— 11^ II < V^d, which completely settles the issue. 
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1. Introduction 

An important problem in the perturbation theory of self-adjoint operators is to study 
the variation of the spectral subspace associated with an isolated spectral subset 
that is subject to a perturbation (see, e.g., [7|). Classical trigonometric estimates in 
subspace perturbation problem have been established by Davis and Kahan [i5 |. For 
further results on subspace variation bounds for self-adjoint operators we refer to 
0, 01, nil, lEl, fl3| and the references therein. 
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In this article we consider a self-adjoint operator A on a separable Hilbert 
space f), assuming that the spectrum of A consists of two disjoint components Gq 
and (7i such that the set ao lies in a finite gap of the set Gi . In other words, we 
suppose that 



where conv denotes the convex hull and overlining means closure. The perturba- 
tions V are assumed to be bounded and off-diagonal with respect to the partition 
spec(A) = ao U di, that is, V should anticommute with the difference E^((7o) — 
E^(o'i) of the spectral projections E^(ao) and E^(<7i) of A associated with the sets 
do and (7i, respectively. For the spectral disposition il.li . it has been proven in |9| 
(see also ifTslfTB l ) that the gaps between Gq and Gi remain open if the off-diagonal 
self-adjoint perturbation V satisfies the (sharp) condition 



where d : = dist(c7oi Ci) stands for the distance between Gq and Gi. Under this con- 
dition the spectrum of the perturbed operator L = A + V consists of two isolated 
components coq C A and coi C M \ A. Here and in the sequel, A denotes the finite 
gap of Gi that contains Gq. (We recall that by a finite gap of a closed set a C M 
one understands an open bounded interval on M that does not intersect this set but 
both ends of which belong to G.) It is worth noting that the norm bound ( 11.21 ) is also 
optimal in the sense that, if it is violated, the spectrum of L in the gap A may be 
completely empty (see IfTOl Example 1.6]). 

The goal of the present paper consists in finalizing a sharp norm estimate on 
the variation of the spectral subspace Ran(E^((7o)) under off-diagonal self-adjoint 
perturbations that was conjectured and partly proven in lfT3l . Our main result is as 
follows. 

Theorem 1. Given a (possibly unbounded) self-adjoint operator A on a separable 
Hilbert space f), assume that its spectrum consists of two disjoint components Gq 
and <7i satisfying condition (II. 11 1. Let V be a bounded self-adjoint operator on Sj 
off-diagonal with respect to the partition spec(A) = Co U <7i and set L — A + V, 
Dom(L) — Dom(A). Assume in addition that V satisfies the bound (11.2b and let 
CO(j ~ spec(L) n A. Then the difference between the spectral projections E^(ao) ond 
El(<Ho) of a and L associated with the respective spectral sets Gq and OJo satisfies 
the norm estimate 



We underline that for 1 1 V 1 1 <d the bound ( 11.31 ) was established in lfT3l . It was 
called there the A priori tan© Theorem. For ||y|| = d this bound may be obtained 
from the result of L13J by continuity. Having proved Theorem[T]we confirm the truth 
of the conjecture of lfT3l Remark 5.7] and we thus close the gap in the subspace per- 
turbation problem for dispositions ( II. Il l which has remained for ||y H/ii € (1, \/2). 
We also remark that the a priori tan Q theorem for eigenvectors 1 3 , Theorem 1.1] is 
a simple corollary of Theorem[T] 



conv((7o) n (7i = and Gq c conv(ai) 



(1.1) 



liVll <\/2c/. 



(1.2) 
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Our proof of Theorem[T]is essentially based on the reduction of the subspace 
perturbation problem under consideration to the study of the operator Riccati equa- 
tion 

XA(^-AiX+XBX = B* (1.4) 

with Ao ^1 ==^|2t[' ^ — ^lai where 2to = Ran(E^((7o)) and 2li — 

Ran(EA(cTi)). In fact, the perturbed spectral subspace £o — Ran(Ei(a)o)) is the 
graph of a particular solution X E ^(2lo,2li) to equation ( 11.4b . In such a case (see, 
e.g., H) 

||EA(cJo)-Ei(a^)|| =sin(ai-ctan||X||). (1.5) 

Thus, having established a bound for the solution X one simultaneously obtains 
an estimate for the norm of the difference of the spectral projections E^(ab) ™d 
Ei(a)o) as well as a bound for the operator angle 



© = ai-ctan\^*X (1.6) 

between the spectral subspaces 2to and £o. For the concept of operator angle and 
related material we refer to [81 and references therein. Note that because of (11.6b 
the operator X itself is usually called the angular operator for the pair of subspaces 
(2lo,£o). 

By (11.5b and (11.6b . the bound (11.3b can be equivalently written in the form 

tan0<M 
a 

which implies that under conditions ( 11.11 ) and ( 11.21 ) the norm of the operator angle 
between 2lo and £o can never exceed the value of arctanV^ (w 54°44'). 

The present article is the third in a series of papers on a priori tan© bounds, 
following li3]|T3|. Its strategy, however, is very different from the approaches used in 
[3 1 and 1 13l. The approach of paper [13] (which was actually the first in the series) 
is based on the properties of sectorial operators and on an involution technique that 
works only in cases where © < 7r/4 (also cf. lO) and the corresponding angular 
operators X in ( 11.5b are contractions. The approach of fi\ only applies to individual 
eigenvectors of L and there is no chance to extend it to multi-dimensional spectral 
subspaces. The key ingredient of the method we use in this paper is a new identity 
for eigenvalues and eigenvectors of the modulus \X \ = \/X*X of X that was found 
only after the articles [31 and [131 were written. Here we mean the identity ( 12.91 ) of 
Lemma l2!2] below which allows us to obtain a norm bound for X even if X is not a 
contraction (see Theorem l3.2l and its proof). 

The paper is organized as follows. In Section|2]we prove the key Lemma IZ21 
Then we recall some known bounds on the shift of the spectrum of the operator A un- 
der a perturbation V satisfying the more detailed (and weaker than (11.2b ) condition 
||y II < -y/c/lAl where |A| stands for the length of the gap A. We also recall a known 
norm bound for the angular operator X in ( 11.5b that is valid for 1 1 V 1 1 < y/d{\A\^^d). 
In Section|3]we employ the identity ( 12. 9b to obtain an estimate for ||X|| already for 
||y|| > Y^^^d^l —d) but in the special case where \X\ is assumed to have an eigen- 
value equal to \\X\\. In Section H] this estimate for \\X\\ is used to prove our most 
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general and detailed subspace variation bound (see Theorem l4.1b . We conclude with 
a proof of Theorem[T] 

The following notations are used throughout the paper By a subspace we al- 
ways understand a closed linear subset of a Hilbert space. The identity operator on a 
subspace (or on the whole Hilbert space) is denoted by Aji^ . If no confusion arises, 
the index 971 may be omitted in this notation. The Banach space of bounded linear 
operators from a Hilbert space 971 to a Hilbert space 91 is denoted by .^(971,91). 
By 971 ® 9t we understand the orthogonal sum of two Hilbert spaces (or orthogonal 
subspaces) 97t and 91. The graph ^(/T) = {y e 97t ® 91 1 y = a; © x e 971} of a 
bounded operator K e ^(971, 91) is called the graph subspace (associated with the 
operator K). By E7-((7) we always denote the spectral projection of a self-adjoint 
operator T associated with a Borel set d C M. The notation p {T) is used for the re- 
solvent set of T . The domain and the range of an operator S are denoted by Dom(5) 
and Ran (5), respectively. 



2. Preliminaries 

It is convenient to represent the operators under consideration as block operator ma- 
trices. Since condition (II. lb will not always be assumed, we first adopt a hypothesis 
that implies no constraints on the mutual position of the spectra of the entries Aq 
and A i . 

Hypothesis 2.1. Let 2lo and 2li be complementary orthogonal subspaces of a sep- 
arable Hilbert space S). Assume that A is a self-adjoint operator on = 2lo ffi 2li 
admitting the block diagonal representation 

^=(1)' m \ Dom(A)=2lo®Dom(Ai), (2.1) 

with Ao a bounded self-adjoint operator on 2to and Ai a possibly unbounded self- 
adjoint operator on 2ti. Suppose that V is an off-diagonal bounded self-adjoint op- 
erator on Sj, i.e., 

V-( 1 f V (2.2) 



B* ^ 

where 7^ B e ^(2li , 2lo), and let L=A + V, Dom(L) = Dom(A), that is, 

L=(g* Dom(L)=2toeDom(Ai). (2.3) 



Under the assumptions of Hypothesis |2T1 an operator X e ^(2to,2li) is said 
to be a solution of the operator Riccati equation (11.4b if 

Ran(X) cDom(Ai) (2.4) 

and (11.4b holds as an operator equality on 2to (cf., e.g., fV, Definition 3.1]). Clearly, 
the solution X, whenever it exists, satisfies X ^0; otherwise, X = implies B — 
which contradicts the hypothesis. In the following by U we denote the partial 
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isometry in the polar decomposition X — U\X\ of X. We adopt the convention that 
U is extended to Ker(X) = Ker(|X|) by 

t/|Ker(X)=0. (2.5) 

In such a case U is uniquely defined on the whole space 2lo (see, e.g., fT, Theorem 
8.1.2]) and 

f/ is an isometry on Ran(|X|) = Ran(X*). (2.6) 

The assertion below provides us with three useful identities for eigenvalues and 
eigenvectors (in case they exist) of the modulus \X\. 

Lemma 2.2. Assume Hypothesis \2.1\ Let X £ ^(2lo,2li) be a solution to the oper- 
ator Riccati equation ( 1 1 A\ . Suppose that \X\ has an eigenvalue X (X > 0) and that 
u, u ^ 0, is an eigenvector of \X\ corresponding to this eigenvalue, i.e. \X\u = Xu. 
IfU is the isometry from the polar representation X — U\X\ of the operator X, then 
Uu G Dom(Ai) and the following three identities hold: 

X{\\AQu\\^ + \\B*uf-\\AiUu\\^-\\BUuf) 

= {l-X^) {{Aqu,BUu) + {B*u,AiUu)), (2.7) 

X{{Aou,BUu) + {B*u,AiUu)) = \\Aou\\^-\\Aouf-\\B*uf, (2.8) 

X^{\\AiUuf + \\BUuf-\\Aouf) = \\Aouf + \\B*u\\^-\\Aouf, (2.9) 

where the entry 

Ao = (/+ |X|2)i/2(Ao+BX)(/+ |X|2)-i/2 (2.10) 
is bounded and self-adjoint on SIq. 

Proof. We start with remark that if X ^0 then Uu = jU\X\u = jXu and, hence, 
Uu G Dom(Ai) by (|231i. For A = we have u G Ker(|X|) = Ker(X) and then Uu = 
G Dom(Ai) by convention ( I2.5l l. We also notice that for the eigenvector m of |X| 
associated with the nonzero eigenvalue A > one automatically has u G Ran(|X|) 
and, thus, in this case the assertion ( 12.6b implies U*Uu = u. 

First we prove the identity ( 12.71 ). If A = 0, then il.li is trivial since Uu — 
due to ( 12.5b . Suppose that A > and set 

( ) = ( U\X\u ) = ( XUu ) ' ^'^■^^^ 

From J/m G Dom(Ai) one concludes that both x and y belong toDom(L). Since X is a 
solution to the operator Riccati equation (11.4b . by, e.g., JT] Lemma 5.3 and Theorem 
5.5] the graphs ^^(X) and ^^(— X*) are reducing subspaces for the operator matrix 
L. Clearly, x G ^(X) and y G ^(-X*) which yields Lx G ^(X) and Ly G ^(-X*). 
Since the subspaces ^^(X) and '^^{—X*) are orthogonal to each other, we have 

{Lx,Ly)^0. (2.13) 
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Using the last equalities in ( 12. lib and ( 12.12b one obtains 

/ Aou + lBUu \ . J ( -IAqu + BUu \ 
^= [ B*u + XAm ) ""'^ \ -XB*u+A,Uu ) ■ ^2.14) 

Substitution of the expressions for Lx and Ly from ( 12.141 ) into the equality ( 12.13b 
results in the identity ( 12.71 ). 

To prove ( 12.8b . we begin with the following equalities: 

AqBX + BAiX =AQBX + B{XA(i + XBX-B*) (2.15) 

^{AQ^BXf-Al-BB\ (2.16) 

by taking into account at the step ( 12.151 ) that, due to ( fTH i, A^X ^ XAq+XBX - B* . 
Since Xu = [/|X|m = XUu and Uu G Dom(Ai), equality (12.16b yields 

X{{Aqu,BUu) + {B*u,AiUu)) = (m, {Ao+BXfu) - ||Aom||^ - ||B*m||^ (2.17) 

Clearly, 

(Ao+BXf = (/+ |X|2)-i/2a2(/+ |X|2)1/2^ (2.18) 

where Aq is the bounded operator given by (12.10b . Since u is an eigenvector of 
by ( 12.18b one obtains 

(m, (Ao+BXfu) = ((/+ |X|2)-i/2„^Ag(/+ |X|2)i/2„) 

= ((1+A2)-1/2„^A2(i+^2^1/2„) 

= (m,A^m). (2.19) 

That the operator Aq is self-adjoint follows, e.g., from fl\ Theorem 5.5]. Hence, 
combining ( 12.17b and ( 12.19b one arrives at ( 12.8b . 

As for the identity (12. 9b . for A = it follows immediately from ( 12.8b . If A > 0, 
then ( 12.91) is obtained by combining (12.8b with (12.71 ). □ 

From now on we assume the spectral disposition ( 11.11) . When necessary, this 
disposition will be described in more detail as follows. 



Hypothesis 2.3. Assume Hypothesis 12.7] Let do = spec(Ao) and Ci = spec(Ai). 
Suppose that an open interval A = (7/, Jr) C M, 7/ < Jr, is a finite gap of the set 0\ 
and Oq C a. Set d = dist((7o, ^i)- 

Below we will use the following assertions obtained by using several results 
proven in ||9l . 



Theorem 2.4. Assume Hvpothesis \2.3\ and suppose that \\V\\ < ^y d\A\. Then: 

(i) The spectrum of the block operator matrix L consists of two disjoint compo- 
nents COo C A and (0\ C K. \ A. /n particular, 

vam{(OQ)>'Yi + {d — rv) and max((Bo) < 7r — (^7— ^v), (2.20) 

where 

(\ 2l|y|| \ 

rv:=||y||tanf -arctan-i^j <t!f. (2.21) 
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(ii) There is a unique solution X G ^(Stoi^li) to the Riccati equation (1 1.41) with 
the properties 

spec(Ao + BX) = 0)0 and spec(Ai -B*X*) = (0\; (2.22) 

the spectral subspaces £o — Ran(Ei(a)o)) and £i = Ran(Ei(oi)) are graph 
subspaces of the form £o — ^(^) <^nd £i = ^^(— X*). 

Remark 2.5. Assertion (i) of Theorem l2.4l follows from 19] Theorem 3.2]. Assertion 
(ii) is obtained by combining |0 Theorem 2.3] with an existence and uniqueness 
resuh for the operator Riccati equation (11.4b estabhshed in Theorem 1 (i)]. 

A sharp a priori norm estimate for the operator angle between the subspaces 
Ran(Eyi((7o)) and Ran(EL(a)o)) and, equivalently, for the corresponding angular 
operator X in (11.5b was obtained in lfT3l Theorem 5.3] under an assumption that 
is stronger than condition ||y|| < -^/t/jAl of Theorem 12.41 We formulate the main 
statement of llT3l Theorem 5.3] in the following form. 

Theorem 2.6 ( II13II ). Assume Hypothesis \2.3\ Assume in addition that 



\\V\\ < ^d{\A\~d). 

Let X be the unique solution to the Riccati equation ( 11.4b with the properties (12.22b . 
Then 

||X|| <tan f^iarctan><(|A|,t/,||y||) j (<l), (2.23) 



where >c{D,d, v) is defined for 
by 

>c{D,d, V 



D>Q, 0<d<^, and < v < y/d{D~~d) (2.24) 



^ ifv<^-^diD-2d), 



2{d{D-d)-v^) 2^ ^ ' 

In the sequel, the estimating function appearing on the right-hand side of 
( 12.231 ) will be denoted by Mi, that is, 

Mi(D,c/,v) :=tan Qarctan><(D,c/,v)^ , (Z),t/, v) e Hi , (2.25) 

where £l\ stands for the set of points (£>,(/, v) e with coordinates D, d, v satis- 
fying ([223. 

Remark 2.7. Using the elementary formula 

(\ \ X 
tan - arctanjc = , x e M, 

V2 ) i + \/T+^' 

one can also write the function Mi (D, li, v) in the algebraic form 

2v 



Ml (Z),t/,v) 



q(0) rf + Vd2+4v2' 



(2.26) 
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Mi{D,d,v) 



where o!!^^ and ClV' denote the corresponding complementary parts of the set Cli 



v{2v+ ^{D-2d)^+4v^) + ^d{D-d) {D-2^d(D-d)) 



Dv+ ^d{D -d)^{D-2d)^+4v^ 



(2.27) 



(1) 



(0) 



(1) 



(ZJ,^/,v)e£2i 
iD,d,v)e£Li 



0< v< '-^d{D^2d)y 
\_^d{p-2d) <v< y'diD-d)^. 



By (12.251 1 we have 

0<Mi{D,d,v) <l for any {D,d,v) eQ.1. 

By representation ( 12.271 ) the function Mi{D,d, v) admits a continuous extension to 
the part 

<9i2i2 |(£>,<i,v) e | £> > 0, 0<d<D/2, v= y/d(p~d)^ (2.28) 

of the boundary of Hi where v ~ y/d{D — d). For the extended function we keep the 
same notation Mi . One verifies by inspection that Mi{D,d,v) = 1 for any {D,d,v) G 
dan- 

Obviously, the function Mi {D,d, v) is infinitely differentiable within the sets 
fl^"'' and Q.[^\ Furthermore, this function and the partial derivatives iMl^M2ll^ 
vary continuously when {D,d, v) passes through the com- 



dMi{D,d.v) 



and 



dM, {D,d.v) 



del 

mon border 



3v 



= ' n of the subsets ilf and Thus, the function 
Ml and its derivatives and are continuous on the whole set £li. 



3. Norm bound for the angular operator in a special case 

Technically, this section is central in the paper We aim at obtaining a norm bound 
for the angular operator X under condition ^t/(|A| —d) < \\V\\ < y^d\A\ which is 
outside of the scope of Theorem l2.6l In the proof we restrict ourselves, however, to 
the special case where the modulus \X\ of X has an eigenvalue coinciding with its 
norm |||X||| = \\X\\. 

In order to formulate the result we introduce another estimating function 

M2{D,d, v) := y 1 + |J -^VdD - ^/{D^d)D - v\ {D,d,v)e^2, (3.1) 
where the set £^2 is defined by 

£22 \{D,d,v)(^9? | £> > 0, 0<d<D/2, y/d{D -d)< v < \/^ | . 

Remark 3.1. Obviously, the function M2 (/),£/, v) is infinitely differentiable inside 
fi2 and continuous on £^2. One verifies by inspection that 

min M2(Z),t/, v) = 1, sup M2{D,d,v) ^ V2, (3.2) 

iD.cLv)£Q.2 (D.d,v)eQ.2 
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and M2{D,d, v) = 1 for any {D,d, v) e dQ.12 where dQ.12 is the intersection (12.28b 
of the boundaries of £2i and Q.2. 



Theorem 3.2. Assume Hvpothesis \2.3\ Assume in addition that 



yJd{\A\-d) < \\V\\ < y/d\A\. (3.3) 

Let X € i^(2lo,2li) be the unique solution to the Riccati equation (ll.4l l with the 
properties (I2.22l l. If\X\ has an eigenvalue jj. such that ji = ||^||, then the following 
bound holds: 

|iX||<M2(|A|,d,||y||), (3.4) 
where the function M2{D, d, v) is given by ( 13. 11 1. 

Proof. Throughout the proof we assume, without loss of generality, that the gap 
A is centered at zero, i.e. Jr ~ — Ji = y\ otherwise, one replaces Aq and A\ by 
Aq = Ao — cl and A\ = A\ — cl, respectively, where c = {yi + jr)/'^ is the center 
of A. The assumption that Ob C A = (—7,7) and d — dist(a(),0'i) > means that 
C7o C [— a,fl] with a = 7— d and ||Ao|| — a. 

Suppose that jU is an eigenvalue of \X\ such that ji ~ \\X\\ = || 1^1 1| and let m, 
||m|| = 1, be an eigenvector of \X\ associated with this eigenvalue, i.e. \X\u — pLu. 
If H = \\X\\ < 1 then, under condition ( 13.3b . the bound ( 13.4b holds automatically by 
the first equality in (13.2b . Further on in the proof we will always assume that jj. > I. 

Let Aq be as in (12.10b . Since spec(Ao) — spec(Ao +BX), from Theorem l2.4l it 
follows that spec(Ao) ~ COq and then (12.20b yields 

0< ||Aom|| <fl + ri/ < 7, (3.5) 
where ry is given by ( 12.21b with |A| = 27 = 2{a + d). At the same time 

\\AiUuf + \\BUuf> \\AiUuf>f, 

taking into account that u E Ran(|X|), ||m|| = 1 and then \\Uu\\ = 1 by ( 12.6b . Hence, 
bvdXSt 

\\AiUu\\^+\\BUu\\^-\\AQuf >f-(a + rvf>0 
and the identity ( 12.91 ) in Lemma IZ21 implies 

^2 ^ \\Aou\\^+\\B*u\\^-\\Aouf 

^ |iAif/M||2 + ||Bf/M||2-||AoM||2- 

Since 

|lAoM|l<fl, |lAif/M|l>7, and < IISlI- 

from ( 13.6b it follows that 

2 . a^ + \\Bf-\\Aou\\^ 
^ - f + \\BUuf-\\Aouf ' 
Because of ;U > 1, from ( 13.7b one infers that 



a 



Bf>f + \\BUuf. (3.8) 
As for the quantity ||Aom||, in view of (13.5b we have two options: either 

0<||AoM||<a (3.9) 
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or 

fl < IIAqmII < fl + rv/. (3.10) 

S — X 

Since for any s, f G M such that t <s the function f{x) : ~ is increasing at < f , 

in the case (|3.9l l combining inequalities (13.7b and (13.8b yields 

In order to treat the case (13.10b properly, one notices that, due to (12.10b . 

||Aom|| =||(/+ i/2(Ao + V2„|| 



<^^-^H||Ao« + MBf/«ll 

= ||AoM||+M||5f/M||, (3.12) 

taking into account that |X|m = pLu at the first step and that \\X\\ = at the sec- 
ond. Since ||Aom|| < a, from ( 13.121 ) one deduces that, in the case (13.10b . ||Bt/M|| > 
^(||Aom|| — a) > and then (13.7b implies 



a^ + \\Bf-\\Kau\\ 
72 + -ij(||AoM||-fl)'-||AoM|l2 



. ::""'r":r'. ... (ifiiAo«ii>«). o.b) 



Inequahty (13.13b transforms into 

||B||2 + 2MK^-|lAo«ll) 
72- ||Aom| 

By combining (13.11b and (13.14b one arrives at the estimate 



M^< """ "JTr 112'' (if||Ao«||>«). (3.14) 



T (p(fl) if IIAomII < fl, 

M <{ Z (3-15) 
^(||Aom||) if ||Aom|| >a, 

where the function <p{z) for z e [0, 7) is defined by 

||B||2 + 2z(fl-z) 
'Pl^)- y2_,2 -■ (3.16) 

One observes that ^(0) = ||B|p/72 > and ^(z) as z — > 7-O since 

||B||2 + 27(fl-7) = ||y||2-rf|A|<0 

by hypothesis (13.3b . Again taking into account ( 13.5b . by (13.15b one concludes that in 
any case 

^2 < max (p{z). (3.17) 

ze[o.y) 

We notice that the function ( 13.16b akeady appeared in the proof of Lemma 2.1 in 
There is a single point zo within the interval [0, 7) (in fact, zo G [0,fl + rv]) where 
the derivative of this function is zero, namely 

if fl = 0, 



^o = <; tMj(fcM)^72 if ,>o. ^'-''^ 
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At zo ths function (p{z) attains its maximum on [0, 7), i.e. 

max (i!)(z) = (i!)(zo). (3.19) 

ze[o.r) 

By inspection, (p{zo) =M2(27,7— a, ||B||)^ ~M2{\A\,d, ||V||)^, where the function 
M2{D,d, v) is given by (13.1b . Combining this with ( I3.17l i and ( 13.191 1 completes the 
proof. □ 

From the two estimating functions introduced in (12.251 1 and ( 13.11 ) we combine 
the total estimating function 



\ Mi{D,d,v) if 0< v< Jd{D-d), 

M(D,d,v):=l ' ' ^ - " (3,20) 

[ M2{D,d,v) if ^d{D-d) < v< Vdl5, 

which is considered on the union H i2i U i22 of the domains i2i and ^2, 

£2= |(D,t/,v) e I D >0, 0<d<D/2, < v < \/^ | . 

Remark 3.3. By Remarks 12.71 and [TTl the estimating function M{D,d, v) is contin- 
uous and uniformly bounded on the whole set £2. It also admits a continuous exten- 
sion to the boundary dO. of Q. (except for the intersection of dQ. with the D axis). 
It should be underlined, however, that the partial derivatives ^^^^^T^' ^'^^gf '^^ 
and ^'^^^•^^'^) are discontinuous when, for d < D/2, the point {D,d, v) crosses the 
common boundary dQ.12 of the sets Cli and £^2. 



4. Subspace variation bound in the general case. 
Proof of Theorem [T] 

The norm bound for the angular operator X obtained in the previous section for 
the special case where \X\ has an eigenvalue equal to ||X|| allows us to prove the 
following general subspace perturbation bound. 

Theorem 4.1. Assume Hvpothesis \2.3\ If\\V\\ < ^/d\A\, then 

||EA((Jo)-Ei((Bo)|| <sin(arctanM(|A|,fl',||y||)), (4.1) 

where (Oq = spec(L) n A and M{D, d, v) is the function defined by (I3.20l l. 

Proof. Assume, without loss of generality, that the gap A lies on the non-negative 
semiaxis, that is, 

< 7/ < 7r; (4.2) 

otherwise, one replaces Aq andAi byA^ =Ao+ (3^ ^c)/ andA'^ =A\ + ^c)I, 
respectively, where c = (7/ + yr)/2. 

First, we consider the case where the spectral subspace 2lo is finite-dimensional. 
Theorem l2.4l (ii) ensures the existence of a unique angular operator X for the pair 
of subspaces 2to = Ran(EA(<7o)) and £0 = Ran(EL((Bo)). Since dim(2lo) < °°, the 
operator X is of finite rank and so is its modulus \X\. Then there is an eigenvalue 
;Uof \X\ such that ^ = |||X||| = ||X||. Hence, for ^/d{\h.\ - d) < \\V\\ < ^^/|A|, the 
bound (EJll follows by (O and (Il20] l from Theorem[32] For ||y || < ^d{\A\ -d) 
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this bound is implied by Theorem l2.6l Therefore, for the case where dim(2to) < °°, 
the bound ( 14.1b has been proven. 

If the subspace 2lo is infinite-dimensional, let {Pi^^jneN be a sequence of 
finite-dimensional orthogonal projections in 2to such that Ran(pi^^) C 2to and 
s-limPi''' = /ain. Using the projections pj"' we introduce the block diagonal op- 
erator matrices 

/ p(0), p(0) „ \ 

A„^{^"^^^" ^ ), Dom(A„):=Dom(A) (= 2loeDom(Ai)), 



which represent the corresponding truncations of the operator A with finite rank 
parts in SIq. We also introduce the finite rank operators 



V„ = 



pi^'B 




and set L„ —A„ + V„, Dom(L„) := Dom(A„) = Dom(A). The operators A„ and Vn 
are self-adjoint. Hence, so are the operators L„. 

Obviously, for any A € C \ M the following operator identities hold: 

(A„-A/)-'-(A-A/)"^ = (A„-A/)"'5„(A-A/)-', (4.3) 

(L„ - A/) - ' - (L - A/) - 1 = (L„ ~lI)-'{S„ + V-V„)iL-lI)-\ (4.4) 

where S„ is the bounded operator on f) given by 

By, e.g., |4, Theorem 2.5.2] we have i-limV,, = V, .s-lim(pi'''AoPi"') = Aq, and then, 
due to (l43b-(l43Tl. 

i-lim(A„-A/)"' = (A-A/)"' and i-lim(L„ - A/)"' = (L- A/)"' (4.6) 

for any A G C\M, which means that both sequences {A„}„gN and {Ln}„en are 
convergent in strong resolvent sense (see, e.g., lfT4l Section V111.7]). 

Let A„ and ¥„ denote the parts of the operators A„ and V„ associated with their 
reducing subspace 

^„=2i|,">®ai, (4.7) 

where 21q"^ = Ran(pi'''). Clearly, the operator A„ is block diagonal with respect to 
the decomposition (14.7b . Dom(A„) = SIq"' ® Dom(Ai), and A„ [^^^ = Ai. Further, for 
the spectral set CTq"' ;= spec (A„U(„)) we have the inclusion 

d^"^ (l[r, + d,yr-d] (4.8) 

and, thus, 

d„ : = dist f spec (A„ | ~ („, ) , spec (A„ | ) j = dist(CT^"' ,ai)>d. (4.9) 
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By its construction, the finite rank operator V„ is off-diagonal with respect to the 
decomposition ( 14.71 ) and 

llK||<||V||. (4.10) 
By the hypothesis we have ||y|| < y^c/]a[. Hence, from (14.10b and ( 14. 9t it follows 
that IIV,,!! < y^c/]a[ < ^/d^^\A\. Then Theorem 12.41 (i) implies that the spectrum of 
L„ : = A„ +V„ consists of two disjoint components Sq""* and such that 

a^o^ c[ri + d,-rl"\Yr-dn + r^yhcA and S{"'cM\A, (4.11) 

where is given by 

Since d < d„ < and V„ satisfies ( 14.101 ), one easily verifies that d„ — ry "^ >d~rv 
with ry given by ( 12.21b . Therefore, from the first inclusion in ( 14.1 lb it follows that 

a),)"' c[Yi + d- rv,Yr-d + rv] for any n e N. (4. 12) 

Furthermore, since the spectral subspace Stg"' — Ran(E^ (Sq"')) is finite- 
dimensional, the bound ( 14.1b applies to the spectral projections (Sg"') and 

||E^^_(CT^"VEi„(S^"')|| <sin(arctanM(|A|,fl'„,||y„||)). (4.13) 

Observing that the function M{D,d,v) is monotonously increasing as the second 
argument decreases and/or the third one increases, by ( |4.9b and ( 14.10b from (14.13b 
one infers that 

||E;j^_(a("VEz„(S("')|| <sin(arctanM(|A|,J,|ly|l)). (4.14) 

Now for an arbitrary e such that < e < li — ry we set := (7/ + e, 7^ — e). Ob- 
viously, by ( 14.8b and ( 14.12b . the open interval contains both sets Sq"' and Wq"'. 

Hence, ^x„(^o^) = E^„(^£) and E^ = ^Z^C^e)- Then inequahty (14041 may 

be rewritten as 

||E^^_(E,)-E£^/E,)|| <sin(arctanM(|A|,«',||y||)). (4.15) 

Clearly, the spectrum of the part L„ I of the operator L„ associated with its 

reducing subspace Jo^ —^jQSjn consists of the single point zero and the same holds 
for the spectrum of the restriction A„ I , i.e. 

spec(L„|_-_^) = spec(A„|-j_) {0}. (4.16) 

By (14.2b this means that none of the sets spcc(L„\~^) and spec(A„|~^) intersects 
the interval E^. Hence, (14.15b yields 

||EA„(E,)-Ei„(E,)|| <sin(arctanM(|A|,^/,||y|l)). (4.17) 
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Meanwhile, equalities (14.16b considered together with the inclusions (14.81) and 
( l4TT2b imply 

{7h7i + d) C p{A„) and {y,- - d,Yr) C p{A„), 

(jhYi + d-rv) C piL„) and (7, - 1/ + rv/, 7,-) C p(L„). 

Then, from the strong resolvent convergence ( I4.6l l of the sequences {A„}„gf} and 
{L„}nsN, it follows (see, e.g., |il4, Theorem V1I1.24]) that for any e such that < 

e <d — rv 

i-limE^„(^£) = Ea(S£) and ,s'-limEz,„(E£) = ElC^e). 

Passing to the limit as « 00 in (|4.17b . one obtains 

||EA(E£)-Ei(E,)|| <sin(arctanM(|A|,J,|lV|l)), 

which is equivalent to ( 14.11 ) since both spectral sets Gq and too are subsets of the 
interval Eg (see Theorem l2.4l (i)) . □ 



Remark 4.2. The bound ( 14.1b is sharp. For ||y || < -\/(i(|A| — d), this has been es- 
tabUshed in IH (see IE] Remark 5.6 (i)]). For y/d{\A\ -d) < \\V\\ < y^d\A\ the 
sharpness of (14.1b is proven by ID Remark 2.3]. For convenience of the reader, 
below we reproduce the corresponding examples from f3\ and fT3l that prove the 
optimality of the bound ( 14.1b . 



Example 4.3 (ID). Let ijo = C and i^i = C^. Assuming that < a < 7 and bi,b2> 
0, we set 

Ao = fl, Ai = ^ ^ ^ ^ , andB= (^1 ^2)- 

The operators (3x3 matrices) A, V , and L are defined on — S)q® = by 
equalities ( 12.1b . (12.2b . and ( 12.3b . respectively. The spectrum Oq = {a} of Aq lies 
in the gap A = (—7, 7) of the spectrum Oi ~ {—7, 7} of A]. Also notice that d = 
dist((705O'i) = 7-fl and |A| = 27. 

First, consider the case where bi—0 and, thus, ||y || = ||Z?|| = b2. In this case, 
for any > satisfying < 27(7+0), i.e., for ||y II < ^|A|(|A| -d) the matrix L 
has a single eigenvalue within the interval A; two other eigenvalues of L are in M \ A. 
For the difference of the eigenprojections Eyi((7o) and Ei{(£)q) we have 



where, as usually, (Oq = spec(L) n A. Since y'lAKlAl — d) > -y/d{\A\ ~ 2d), equal- 
ity ( 14.18b proves the optimality of the bound ( 14.1b for 

||y||<lv/^(|A|-2of) 

(see definition ( 12.26b of the restriction A/ |^(o) = Mi |^(o) of the function M onto the 

set nf'). 



Second, assume that bis a positive number such that 



y/y2-a2<b<^2y{y-a) (4.19) 
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and let (cf. formula (13.181 1) 

if fl = 0, 



Zo:^{ 27^-^2 f2f-b^\' , ^ (4-20) 

■ ' ' ' - y2 if a > 0. 



2a V \ 2a 

The second inequality in (|4.19t implies zo E A. Note that, under condition ( 14.191 ). 
for 

one has < f < 1 (actually, f < 1/2). Then set 

bi:=^/Y~tb, b2:^Vtb. (4.21) 

If bi and b2 are introduced by (14.21b . then ||y j| =b. Furthermore, in this case ( |4.19t 
is equivalent to 

.yd{\A\-d) < \\V\\ < ^/d\A\ (4.22) 
and the number zo given by (I4.20l i represents the single eigenvalue of the matrix L 
within the interval A, that is, (Oq = spec(L) n A = {zo}- An explicit computation of 
the only eigenvector of L corresponding to the eigenvalue zo results in 

||EA(c7o)-EL(a)o)|| =sin(arctanM2(|A|,t/,||y||)). (4.23) 

Taking into account definition ( I3.20l i of the function M, equality (14.231 1 proves the 
sharpness of the bound ( 14. Il l for V satisfying (I4.221 I. 

Example 4.4 ( lfT3l ). Leti3o=-^i = C^. Assume that Ao,Ai, andZ? are 2 x 2 matrices 
given respectively by 

where < a < 7 and bi,b2 > 0. Let the operators (4x4 matrices) A, V, and L 
be defined on = ^jo(BSji = by equalities (12. 11 1. ( 12.2b . and ( 12.3b . respectively. 
Clearly, the spectrum Gq = {—a, a} of Aq lies in the gap A = (—7, 7) of the spectrum 
C7i {—7,7} of Ai. We also note that d — dist((7o,0'i) = J- a, |A| = 27, and ||y|| = 

\\B\\=bi+b2. 

One verifies by inspection that the 2x2 matrix 



X=[ ^' ^2 

-H2 -Hi 



where 



Hi = 



K2 



2b\^{y+af+Ab 



{y+a)^{y-af + Ab\ + {y-a)^{y+af + Abl 



{y+a)^{y-af+Ab\ + {y-a)^{y+af+Abl 

is a solution to the operator Riccati equation ( 11.4b . Moreover, under condition ||y || < 
^/t/lAj = \/2y{y— a) the spectrum of Aq + BX lies in the interval A, while both 
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eigenvalues of Ai —B*X* are in M\A. Applying, e.g., |f9, Theorem 2.3] one con- 
cludes that the graph subspace ^(X) = {x©Xx| x G Sjo} is the spectral subspace of 
L associated with the spectral set cOq = spec(L) n A. By (11.5b this yields 

||EA((7o)-EL(a\))|| =sin(arctan(>fi +><2)), (4.24) 

taking into account that \\X\\ = ki+ >^2■ 
Now pick up an arbitrary b satisfying 

^^2{y~a)a<b<^/y^-a^ (4.25) 



and set 



where 



1 /, . , 1 



p) eindb2^-{b-p), (4.26) 



2\ ^ 2' 

if fl = 0, 



a 



- [^Y^t>^ + a^{f-a^-b^)-rb) if a > 0. 



The positivity of bi is obvious. The positivity of b2 is implied by the first inequality 
in ( 14.251 1. Since ||y|| = = fci + fc? = fe. the pair of inequalities (14.25b is equivalent 
to 

^-^d{\A\-2d)< \\V\\ < ^d{\A\-d). (4.27) 
Evaluation of the sum k\ + xi for b\ and b^ given by ( 14.261 ) reduces ( 14.241 ) to 

||E^((Jo)-El(cOo)|| -sin(arctanM^(i)(|A|,d,||y||)), (4.28) 

where M (1) :=Mi| (i) =A^L(i) is the restriction ( l2.27b of the function M onto the 



'I 



set i2j . Therefore, equality ( 14.281 ) proves the sharpness of the bound ( 14.1b for V 
satisfying ( 14.27b . 

Theorem[T|is nothing but a corollary of Theorem l4. 1 1 

Proof of Theorem^ Set 2to = Ran(EA(c7o)) and 2li = Ran(EA(cTi)). With respect 
to the orthogonal decomposition = 2lo 2li the operators A and V are block 
operator matrices of the form ( 12.11 ) and (12.2b . respectively. The length of the gap 
A satisfies the estimate |A| >2d and, hence, condition ( 11.2b implies ||y|| < ^(i|A|. 
Then by Theorem l4.1l we have estimate (14.1b . It remains to observe that, given the 
values of ||y|| and d satisfying ( 11.2b . M{D,d, \ \V\\) is a non-increasing function of 
the variable D, D> 2d. For D varying in the interval [2d,°o) it attains its maximal 
value aiD~2d and this value equals 

max M(D,d,\\y\\)^M{2d,d,\\V\\^ 

D:D>2d d 

Hence, (14.11 ) yields ( 11.3b . completing the proof. □ 



Remark 4.5. Example 2.4 in [|3j (representing a version of Example 14. 3l for a = 
and b\=b2 = b/^/2, where b>0) shows that the bound (11.31 ) is sharp. 
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